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Abstract 

In light of recent findings which seem to disfavor a scenario with (warm) dark matter entirely 
constituted of sterile neutrinos produced via the Dodelson-Widrow (DW) mechanism, we investi- 
gate the constraints attainable for this mechanism by relaxing the usual hypothesis that the relic 
neutrino abundance must necessarily account for all of the dark matter. We first study how to 
reinterpret the limits attainable from X-ray non-detection and Lyman-a forest measurements in 
the case that sterile neutrinos constitute only a fraction f s of the total amount of dark matter. 
Then, assuming that sterile neutrinos are generated in the early universe solely through the DW 
mechanism, we show how the X-ray and Lyman-a results jointly constrain the mass-mixing pa- 
rameters governing their production. Furthermore, we show how the same data allow us to set a 
robust upper limit f s ^ 0.7 at the 2a level, rejecting the case of dominant dark matter (f s = 1) at 
the ~ 3a level. 

PACS numbers: 95.35.+d,14.60.Pq,14.60.St 
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I. INTRODUCTION 



Recent astro 
galaxies [l], 0, 
sized galaxies 



sical observations, such as the indications of central cores in low-mass 



0, i| , and the low number of satellites observed in Milky Way- 
indicate two possible shortcomings of the [A]CDM paradigm 



see |l3l . Il4| ) and have boosted the interest in a warm dark matter (WDM) scenario which 




may alleviate these possible small scale problems [l5|, Electroweak singlet right handed 
( "sterile" ) neutrinos v s with mass in the keV range appear to be appealing candidates since 



18, QJ, 20, 21, 22 of the Standard Model 



they naturally arise in many extensions [Tj 
and they could be produced in the early universe through the Dodelson-Widrow (DW) 
mechanism [23[ involving (non-resonant) oscillations with the active species. 

Direct constraints on mass and mixing ang le of the sterile neutrino 1 can be obtained by 



exploiting X-ray observations [261 . 1271 . l28l . l29j |. Indeed, sterile neutrinos posses a radiative 



decay channel |30l . l3l| , which gives rise to photons potentially detectable in astrophysical 



X-ray and 7-ray sources. Limits of such kind have been derived exploitin g di ffuse cosmic X 



ray background [27. l32l. M ilky Way "blank sky" observations [33|, ]34], ]35|, |36j , X-ray spectra 



38l . l39l . |40[ , and gamma-ray line emission 



of nearby galaxies [35|, ]37J or clusters of galaxies 
searches in the galactic center region [41]. Indirect constraints on the neutrino mass can 
be achieved by the observations of cosmological structure formation on very small scales, 
where WDM typically suppresses the clustering process. In this context, Lyman-a forest 
measurements (henceforth Ly-a) constitute the most suitable probe (42|, being sensitive 
tracers of the primordial density fluctuations on the smallest scales. 

Assuming that sterile neutrinos are produced solely through the DW mechanism and 
requiring that they must account for all the dark matter density Q s = fW, the recent X-ray 



analyses [32|, ]33|, [34 



371 . l38l . l39l . |40j determine upper limits on their mass which lie 



in the range (3-8 keV). On the other hand, under the same hypotheses, the latest Ly-a 



analyses [43j, |44| of the high redshift flux power spectra, measured by the SDSS survey [45 



furnish lower bounds in the range (10-13 keV), in clear tension with the X-ray upper limits. 2 
While these results imply that sterile neutrinos cannot account for the entire amount of 
dark matter in the universe, 3 the possibility exists that the DW mechanism could have pro- 
duced only a fraction f s = Q s /Q dm of its total content. This does not prevent sterile neutrinos 
from playing an important role in other contexts, such as pulsar kicks 



supernova explosions [57], 58| 




56], 
65]. 



and reionization of the universe [59l . 
Indeed, it has been already recognized [49| that a sterile neutrino explanation of the pulsar 
kicks can be reconciled with the existing X-ray constraints assuming the (lowest possible) 
relic abundance provided by the DW production mechanism. Furthermore, the suppression 
of small scale structures induced by WDM, while reducing problems at low redshifts, de- 



lays the onset of reionization (59I . 60| . possibly to an extent difficult to reconcile with the 

Whether the X-rays emitted by sterile neutrinos can effectively di- 

and a sub dominant (f s < 1) 



minish this inconsistency is still an open question [62j, ]6J, [65 



WMAP3 findings [66 



24j and MiniBooNE 25j experiments are not sensitive to the extremely small 



1 Notice that the LSND 
mixing angles involved in the DW mechanism. 

2 Similar lower bounds (> 10 keV) have been determined by QSO gravitational lensing observations 461 ] . 

3 This conclusion can be evaded invoking alternative production mechanisms (not considered in this work) , 



such as resonant oscillations 



reheating temperatures 5(J, ]51 1 



471 ] . generation operative during [48| or prior to 49( oscillations or low 
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relic abundance would probably reduce the tension also in this context. 

Therefore, it seems interesting to investigate in a systematic way the constraints attain- 
able on the mass-mixing parameters which regulate the DW production mechanism relaxing 
the usual assumption of "dominant" dark matter (fl s = O dm ). With this purpose in mind, 
we first address the issue of reinterpreting the X-ray and Ly-a constraints in the presence 
of a u s subdominant relic abundance (Q s < fidm) which is fixed a priori, i.e. independently 
of the oscillation parameters. Then, we treat the general case in which the relic abundance 
depends on mass-mixing parameters as predicted by the DW mechanism. We determine 
the range of such parameters that is compatible with the restrictions imposed by the joint 
analysis of X-ray and Ly-a measurements. In addition, we show how these results indirectly 
compel the relic abundance to be smaller than f s < 70% (at the 2a level), even when the 
large (hadronic) uncertainties affecting the theoretical calculations are taken into account. 

The paper is structured as follows. In Sec. II, we present our constraints obtained from the 
diffuse X-ray background, briefly discussing their behavior for values of f s < 1. In Sec. Ill, 
we introduce a recipe to derive lower bounds on the sterile neutrino mass for values f s < 1, 
appropriately rescaling the Ly-a limits already existing for the usual case f s — 1. In Sec. IV, 
we show how the tension existing between X-ray and Ly-a constraints is reduced when the 
sterile neutrinos are subdominant. In Sec. V, we discuss the restrictions imposed on the 
mass-mixing parameters by the joint analysis of X-ray and Ly-a data, also determining the 
maximum value of the relic abundance allowed by current data. Finally, we present our 
conclusions in Sec. VI. 



II. CONSTRAINTS FROM DIFFUSE X-RAY DATA 

Sterile neutrinos possess a (subdominant) radiative decay channel v s — > ^v a into an active 
neutrino v a and a photon 7 with energy equal to a half of the sterile neutrino mass E = m s /2. 
For Majorana neutrinos 4 the radiative decay width can be expressed as 5 

where Gp is the Fermi constant, a is the fine-structure constant and 9 is the mixing angle 
in vacuum between the sterile neutrino and the active species which, for the very small 
values considered in this work, is related to the effective mixing angles 9 a 's with each of the 
standard (active) neutrinos as, 

sin 2 6 ~ 6 2 = €■ (2) 

a=e,fx,r 

In order to derive constraints from the diffuse X-ray background we use the data provided 
by the HEAO-1 detector [67|,|68| since this experiment covers a large energy band, allowing us 
to put constraints up to large values of the mass (~ 50 keV). This becomes relevant when the 
hypothesis that sterile neutrinos are dominant is relaxed. Furthermore, these measurements 
have been shown to furnish stringent and robust constraints [32| , which are comparable with 



those obtained from nearby galaxies [35|, |37|] or clusters of galaxies [38|, |39|, |40 
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For a Dirac sterile neutrino the decay width is a half 31] 



5 Throughout the paper we use natural units with Ti = c = ks 
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The signal expected in the HEAO-1 detector due to the decay of the putative dark mat- 
ter candidate originates in part from the extra galactic (EG) neutrinos, and in part from 
neutrinos clustered in the Milky Way (MW) halo. The EG contribution can be evaluated 
assuming a uniform distribution of neutrinos in the visible universe up to very small red- 
shifts. Defining Fe as the present energy flux of photons produced by neutrino decays, the 
differential energy flux (energy flux per unit energy and solid angle) can be expressed as [69| 

J2 77-EG p o _ 

EG _ " E _ r 1 7 "dm Pc /g\ 

Ve ~ dVLdE Js A^m s H{m s /2E -I) ' {) 



where p c is the present critical density, fldm — 0.20 is its fraction in form of dark matter [66 
and H(z) is the Hubble function which, assuming a flat A-matter dominated universe, is 
related to the present expansion rate 6 H as 



H(z)~H ^n A + n m (i + z)3, (4) 

where fl m ~ 0.24 and — 0.76 are the present fractions of the critical density [66] in form 
of (total) matter and dark energy respectively. 7 

The monochromatic energy flux produced by the sterile neutrinos clustered in our galactic 
halo can be formally expressed as 

,.mw _ <^£e_ _ , T 7 m s f m s \ 

Ve = imE - fs ^2E Sdm 6 r " ~TJ ' (5) 

which depends on the direction through the mass column density Sdm of dark matter along 
the line of sight (l.o.s.), 

Sdm = J pdm(x)dx . (6) 
l.o.s. 

It is useful to compare the relative magnitude of the EG and MW contributions through the 
(directionally-dependent) ratio of their fluxes integrated over energy. From Eqs. (EH) one 

gets 

Ul °- S - j^dE ~ ^10- 2 g cm" 2 J ' 1 ' 

which, for typical mass column densities Sdm, turns out to be of order unity. In particular, 
adopting the recent evaluation of the galactic halo mass determined in 170] for a Navarro- 



Frenk- White profile, we find values of R in the 2a range (0.5-3.0), in the direction of the 
galactic anti-center (corresponding to the lowest column density). In consideration of the 
all-sky coverage of the HEAO-1 detector 67, 68[, we adopt the simple and conservative 



choice R — 1, thus assuming a MW contribution independent of direction. The effect of 
a more realistic treatment (requiring precise knowledge of the sky coverage of the detector 
during its time of operation) would only render our limits more stringent. 

In order to find constraints on the mass-mixing parameters, we have performed a spectral 



analysis of the HEAO-1 data following the procedure described in [32|, adding the EG and 



6 H ~ 100 h km s" 1 MpcT 1 with h = 0.73 [66 

7 We have checked that the uncertainties on the cosmological parameters have a negligible impact on our 
results. 
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MW contributions to the continuous X-ray spectrum parameterized by the empirical formula 



provided in [68|. In Fig. 1 we show the region excluded at the 3a level (above the curves) 
obtained in the case f s = 1. For completeness, we show separately the constraints that 
we obtain including only the EG contribution (dotted line), and only the MW contribution 
(dashed line), together with the constraints obtained properly taking both contributions into 
account (solid line). Even if the total MW and EG fluxes have the same order of magnitude 
(exactly identical in our case), the MW contribution has a sensibly larger impact in deter- 
mining the constraints. This different sensitivity is due to the different form of the energy 
spectrum produced in the two cases. Indeed, while the MW signal is just a line broad- 
ened into a (large) gaussian (centered around E = m s /2) by the (poor) energy resolution 
(AE/E ~ 25%) of the detector j68|, the EG spectrum is intrinsically broadened towards 
low energies as a result of the integration over redshift and is further enlarged by the effect 
of the resolution of the detector. Both constraints obtained with or without the inclusion 
of the MW contribution are in good agreement with the analogous ones found respectively 
in [32| and [2^] . Our bounds appear only slightly weaker when the MW contribution is also 
included, presumably due to the different choice adopted for the mass of the galactic halo. 

As is evident from Eqs. (JIHSD, the differential flux produced by the decay of sterile neu- 
trinos depends on the product f s sin 2 29. Hence, for a fixed neutrino mass, the upper limit 
on the mixing angle must weaken for decreasing values of f s , rescaling exactly as the inverse 
of the fraction f s . This means that if one fixes a priori the fraction f s < 1 (i.e. indepen- 
dently of the mass-mixing parameters), the excluded region in the plane [m s ,sin 2 2#] will 
reproduce that found for the usual case f s = l modulo a "rigid" shift towards larger values 
of the mixing angle (see Sec. IV). The form of the contours will be altered in a different way 
when, as predicted by the DW mechanism, the relic abundance depends on the mass-mixing 
parameters (see Sec. V). 



III. CONSTRAINTS FROM THE LYMAN-a FOREST 

Massive neutrinos affect the cosmological evolution through two distinct effects. Firstly, 
in the early universe, they are relativistic and therefore contribute to the overall radiation 
budget of the universe, shifting the epoch of matter-radiation equality. Secondly, neutrinos 
erase fluctuations on scales smaller than the horizon size at the epoch when their kinematics 
becomes non-relativistic - the so-called free-streaming length [7l|. For neutrinos under dis- 
cussion here, the first effect is completely negligible while the second one can leave important 
imprints on structure formation which are detectable by Ly-a forest measurements. Indeed, 
these observations are directly sensitive to the free-streaming length which can be roughly 



expressed as [72 



where (p/T) is the mean momentum over temperature for the sterile neutrino distribution. 
To impose constraints at such low length scales (< Mpc), one needs a probe of the small-scale 
matter power spectrum at high redshift, where the information on the primordial fluctuations 
has not yet been completely lost due to non-linear evolution. The Ly-a forest at redshifts 
2-4 satisfies both requirements, constituting the most suitable probe of WDM j42|. Two 
recent and independent Ly-a analyses have provided very stringent lower bounds on the 
neutrino mass. Assuming that all the dark matter is made of sterile neutrinos and that they 
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are produced via non-resonant oscillations, a lower bound of 13 keV (at 95% C.L.) has been 
determined in (43)], and a slightly weaker limit of 10 keV (at 95% C.L.) has been obtained 
in 0. 

The most appropriate approach to constrain mixed models with f s < 1 (we assume that 
the remaining fraction (1— f s ) of dark matter is made of CDM) would be by running a grid of 
hydrodynamical simulations in order to relate the observed Ly-a flux power spectrum with 
the expanded set of parameters considered here. This exceeds the scope of this work, and 
we instead adopt a different method which allows us to obtain reliable constraints simply 
rescaling the ones found in [43j for the case f s = 1. 

Assuming fiducial values for all cosmological parameters, we calculate the growth of 
perturbations including a mixed contribution of WDM (/ s ) and of CDM (1 — f s ) in the code 
CAMB [73[], incorporating the v s relic abundance through the relation 

Vsh 2 = (3 (-^-) , (9) 



93.2eV 

where (3 is a suppression factor and we assume that sterile neutrinos posses a thermal 
momentum distribution. Recent calculations 74], |75j have shown that the real distribution 
exhibits appreciable deviations from the thermal form, which can be rou ghly approximated 
by a moderate shift of the average momentum toward lower values 74|, [75|. This implies 
[see Eq. flH])] smaller free-streaming lengths and a consequent weakening of the lower bounds 
on the mass derived with the assumed approximation of a thermal distribution. The limit 



of 13 keV quoted in [43|] already accounts for a 10% correction {(p/T) ~ 0.9) as evaluated 



in 



741 ] . Since the calculations [751] indicate deviations as large as ~ 20% ((p/T) ~ 0.8 for 
masses of 10 keV), we have conservatively rescaled the lower bound obtained in [43| by a 
further factor 0.9, thus adopting for the pure WDM case (/ s = 1) the reference lower bound 
m s > 11.5 keV. 

Due to their pencil-beam nature, Ly-a measurements effectively probe the projected ID 
flux power spectrum rather than the full 3D power spectrum. This implies that small-scale 
modes are projected onto much larger modes (43[, with enhanced sensitivity on the scales 
of interest. Therefore, for each model (defined by the values of f s and m s ), we integrate the 
resulting 3D linear power spectrum into the relevant ID spectrum, using the relation 

1 r°° 

Piv(k) = — / P 3D (k)kdk . (10) 

ZTT Jk 



In order to derive lower limits on the neutrino mass from those already calculated in [43 
for the case of pure WDM (f s = 1), we adopt the following rescaling procedure (76[. For a 
model with a WDM fraction f s < 1, we find the value of the neutrino mass that produces 
a suppression of the ID power spectrum at the fiducial pivot wavenumber kf = 2h Mpc -1 
equal to that one obtained (for a larger value of the mass) in the case f s — 1. Once the 
confidence level (C.L) at which a given mass is disfavored in the pure WDM case is provided, 
we assume that the values of the mass calculated as above for the models with f s < 1 are 
disfavored at the same C.L. 

Our procedure can be visualized in Fig. [2] which displays in the plane [m s , f s ] the iso- 
contours of the fractional suppression of the ID matter power spectrum obtained for mixed 
(CDM +WDM) models with respect to the case of pure CDM. The values refer to the pivot 
scale k — 2h Mpc -1 . The region under the thick solid line is excluded at the 2a level 
and corresponds to a fractional suppression of ~ 23%, which is the value found in [43| for 
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the case of pure WDM corresponding to the quoted 2a lower bound on the neutrino mass. 
We observe that for small values of the relic abundance the lower limits depart from the 
power-law behavior displayed at large values, indicating a loss of sensitivity of the Ly-a 
measurements. 

Although we have chosen the pivot wavenumber kf = 2h Mpc -1 as this is the scale where 
SDSS Ly-a data are most sensitive |j43j], we have checked that our rescaling procedure returns 
stable results with respect to possible different choices of the pivot wavenumber in the range 
(1-5 h Mpc -1 ). In fact, a different choice of kf in such range alters the lower limits by less 
than 10%, provided that sufficiently high values of the WDM fraction (f s > 0.1) are consid- 
ered. For smaller values of f s , the discrepancies become more pronounced and our procedure 
should be replaced by a more accurate treatment involving detailed simulations. We also 
observe that our rescaling procedure tacitly assumes that the change of cold to warm dark 
matter fraction affects just the linear power spectrum while leaving the transfer functions 
from the linear ID power spectrum to the observed flux power spectrum unchanged. Cor- 
rections to this simple approximation could arise since the ID power spectrum probes very 
high k modes which are probably significantly affected by the Jeans smoothing, non-linear 
coupling from larger scales and redshift-space distortions [77|. Therefore, dedicated hydro- 
dynamical simulations are required in order to render our conclusions more quantitative. 



IV. THEORETICAL PREDICTIONS AND EXPERIMENTAL CONSTRAINTS 

In this section we will show how the constraints coming from X-ray and Ly-a data relate 
to the theoretical restrictions imposed by non-resonant production by confronting the usual 
case of dominant dark matter (f s = 1) with a representative case in which it is subdominant 
(f s = 0.2). The results of this sections will serve also as a guide for the interpretation of the 
more general results presented in the next section. 

For our purposes we use the theoretical calculations performed in [75j. Although similar 
calculations have been recently performed j74[, which are in substantial agreement with those 



presented in [75], the latter incorporate the most conservative evaluation of the hadronic 
uncertainties affecting the production process. Furthermore, fitting formulae valid in the 
whole range of mass and mixing parameters explored here are provided which allow us draw 
iso-abundance curves in the mass-mixing plane and to incorporate the theoretical predictions 
in the analysis. According to the calculations performed in [75[ , after a slight change in the 
notation, the relic abundance of sterile neutrinos can be conveniently expressed as 

^owi^^vsvfm (id 



where the first term in parentheses embeds a mild dependence on the neutrino mass and on 
the flavor structure 8 through the slowly varying functions C a , which are directly evaluated 
in [75j]. From direct comparison of Eqs. (Q and (TIT]) , we observe that the suppression factor 
/3 in Eq. (Q is related to the production mechanism. In particular, once the relic abundance 
and the mass are fixed, the factor (3 is unequivocally determined by the mixing angles. This 
dependence, together with the deviations from the thermal distribution, distinguishes the 
case of (out-of-equilibrium) production via non-resonant oscillations from the case of early 



The functions C a exhibit a moderate hierarchy [75| (C e > C M > C T ). 
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decoupling (in equilibrium), where an analogous suppression factor is instead related to the 



number of relativistic degrees of freedom at the time of decoupling [78 



The fitting formulae presented in [75] provide the relation that mass and mixing param- 
eters must satisfy if sterile neutrinos constitute all the amount of dark matter. These are 
derived imposing Q s = f2 dm in Eq. (fTTl) . and thus correspond to the case f s = 1. Given the 
form of Eq. fill I) , one can obtain fitting formulae for the fractional relic abundance simply 
rescaling those ones valid for the case f s = 1. In particular, for the best estimate of the 



hadronic effects, one gets the average abundance [75 



logio (fD 



( m s \ , / sin 2 20' 
+0.17+1.84 log 10 y,)+log 10 ^ 



(12) 



which we will use as the reference estimate in our analysis. Taking into account the hadronic 
uncertainties, the authors of 75j determine two extreme cases, corresponding respectively 
to the minimal 



logio (fi 



-0.07 + 1.74 lot 



10 



m s \ 1 

keV ) + l0gl ° 



'sm*26 s 

io- 7 . 



(13) 



and to the maximal 



logio (fs 



( m s \ , ( sin 2 20' 
+0.62 + 1.74 logio (^) +logi 



410 



V io- 



(14) 



abundances that can be produced for given mass-mixing parameters. These estimates consti- 
tute two extreme cases in that they have been obtained by "pushing" all the errors induced 
by different and in principle independent sources of uncertainty 10 towards the same direction. 
Lacking a precise definition of their statistical significance, we will adopt the conservative 
choice to interpret them as 2a limits in our analysis. 

In order to show how the constraints imposed by X-ray and Ly-a data relate to the theo- 
retical predictions, we superimpose in Fig. [3] all the (2a) experimental restrictions obtained 
for the two fixed values f s = 1 (left panel) and f s = 0.2 (right panel) to the theoretical curves 
corresponding to these same values and obtained by Eqs. (11211141) . In both panels, the region 
above the descending thick solid curve is excluded by the X-ray analysis while the region 
below the thick horizontal line is forbidden by Ly-a. The inclined dashed line corresponds 
to the best estimate given by Eq. (Jl2"]h while the other two dotted lines represent the two 
extreme cases of Eqs. (j!3til4p . 

In the left panel, corresponding to the usual case / s = 1 (see Refs. [37J, |44j for analogous 
plots), the strong tension existing between experimental constraints and theoretical predic- 
tions is clearly evidenced: the theoretical curves lie in the region which is excluded by X-ray 
data, by Ly-a data or by both sets of data. In particular, looking at the intersections points 
between the theoretical curve and the border of the regions excluded by X-ray and Ly-a 
data one find the incompatible limits m s < 7 keV and m s > 11.5 keV respectively. Hence, 
as expected, the dominant case f s = 1 seems strongly disfavored. 



Here, we have taken the average of the two cases respectively dubbed as "case 1, mean" and "case 2, 
mean" in 



75J, obtained for two extreme flavor structures. 



10 These are mainly the uncertainties on the equation-of-state (EOS) and on the hadronic scattering pro- 
cesses [75i | . 
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In the right panel, all the curves (including the theoretical one) refer to fractional abun- 
dance f s = 0.2. Here the iso-abundance lines move (for a fixed mass) toward noticeably 
smaller values of the mixing angle. On the contrary, as discussed in Sec. II, the X-ray limits 
become less stringent "moving" toward larger (by a factor l/f s = 5) values of the mixing 
angle, while the Ly-a lower bound on the mass drastically decreases (m s ^2.5 keV) as 
discussed in Sec. III. The net result is that for subdominant dark matter with fractional 
abundance f s = 0.2 a region for the mass-mixing parameters is permitted. If one adopts the 
best estimate for the relic abundance, this region corresponds to the segment of the dashed 
line delimited by the two regions excluded respectively by X-ray and Ly-a data. The end 
points of such segment correspond to values of the neutrino mass in the range (2.5-16 keV) 
and to values of the mixing angle in the range (9 x 10 -11 -2.5 x 10 -9 ). Intuitively one can 
expect that such ranges would change (reducing) for increasing values of the fractional abun- 
dance f s until they disappear for sufficiently high values of f s . This is the case as we will 
show by the quantitative analysis performed in the next section. 



V. CONSTRAINTS ON THE DODELSON-WIDROW SCENARIO 

In this section we will show how the results coming from X-ray and Ly-a data can be 
used to put constraints on the parameters which govern non-resonant production of sterile 
neutrinos in the early universe. In order to determine the mass-mixing parameters which are 
compatible with the experimental data, we impose that for each point of the mass-mixing 
plane, the relic abundance is the one determined by the theoretical calculations. 

Figure 0] shows the results of such procedure for the three different estimates provided by 
Eqs. fll2H14p . The first panel corresponds to the "average" case [Eq. (fT2l) ]. while the second 
and third panels represent respectively the two extreme cases of lowest [Eq. (fl3|) ] and highest 
[Eq. ffT4j) ] production efficiencies. In each of the three panels the solid curve delimits the 
parameter region (below) allowed at 2a level, while the dotted lines represent iso-abundance 
contours. The left upper branch of the curve is determined by the X-ray constraints which 
now display a different behavior with respect to those presented in Fig. [3J Indeed, the relic 
abundance is not fixed a priori but varies with the parameters. Consequently, the limits 
on the neutrino mass now decrease faster with the mixing angle due to the corresponding 
increase of the relic abundance along the curve. As already observed in |49[ these limits 
are less stringent than those obtained for the case f s = 1 (see the left panel of Fig. [3]), 
thus leaving room for a sterile neutrino explanation of the pulsar kicks. The allowed region 
is limited on the right by the fast descending branch determined by Ly-a measurements 
which, fixing a lower bound on the mass, indirectly limits also the mixing angle. The two 
branches converge for increasing values of f s "closing" the allowed region in correspondence 
of a maximum value of the relic abundance (evidenced by a dashed segment) which depends 
on the assumed theoretical calculation. We observe that the left upper branch of the allowed 
region in the top panel of Fig. H] is similar to that one determined in |49| using the same 
procedure. However, since only the X-ray data are considered in the right lower branch 



of the allowed region is determined in [49J] by the "overclosure" condition (f s = 1), which is 
less restrictive than the constraints we obtain including also the Ly-a data. Moreover, we 
stress that only the inclusion of both set of data (X-ray and Ly-a) is able to jointly constrain 
the fractional abundance. The comparison of the three panels of Fig. H] allows a qualitative 
evaluation of the impact of the theoretical uncertainties. In particular, the allowed region is 
subject to appreciable changes, enlarging (reducing) for lower (higher) production efficiencies 
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and the upper limit on the abundance is subject to noticeable fluctuations in the range (0.55- 
0.75). 

In order to incorporate such uncertainties in a quantitative way in our analysis we al- 
low the fractional abundance f s to vary around its average value [determined by Eq. (TT21 ] 
according to a log-normal 11 distribution with standard deviation taken equal to one half of 
the excursion (respect to the average) determined by the two extreme cases in Eqs. (I13H14I) . 
This corresponds to adding a penalty factor 12 rj to the total x 2 

X 2 = Xx-my( sin2 26 '> /«) + xt y - a (™s, fs) + m s, fs) , (15) 

defined as 

r io gl0 (/ s )-iog 10 (/r) 

' L Alog 10 (/ fl ) 
with la (asymmetric) errors given by 

A \og 10 (f s ) = 0.5 [log 10 (/r x ) - logio(/D] (/. > fT) , (17) 

A log M (/.) = 0.5 [log 10 (/r) - log 10 (/r n )] (/* < fD • (18) 

Then, in order to obtain the allowed region for the mass-mixing parameters, we marginalize 
the x 2 m Eq. f|T5|) with respect to the parameter f s . Figure [5] shows the allowed region 
(at the 2a and 3a levels) obtained following this procedure. As expected, the 2a region is 
slightly enlarged respect to that one obtained for the "average" case (first panel of Fig. H]). 
Note that in this plot iso-abundance contours cannot be drawn since the value of f s is now 
determined by the marginalization process and thus depends on the confidence level. In 
particular, the maximum value allowed at the 2a level turns out to be f s < 0.7 which is 
slightly more stringent than the absolute upper limit (0.75) obtained for the extreme case of 
the highest production efficiency (third panel of Fig. H]). In fact, any departure from the best 
estimated abundance is now hindered by the counterbalancing effect of the penalty factor 
rj. Furthermore, we find that the case in which sterile neutrinos constitute all the amount 
of dark matter (f s = 1) is disfavored at the ~ 3a level. 

As already noted in the previous section (see the discussion concerning the case f s = 
0.2), for a given value of the relative abundance f s produced via non-resonant oscillations, 
the experimental data determine two allowed ranges for the mass and the mixing angle 
respectively. From Fig. H] one can see that the amplitude of such ranges gradually diminishes 
until it eventually becomes zero for the maximum value allowed for the abundance. This 
behavior can be effectively visualized marginalizing the x 2 m Eq. ffT5l) with respect to one of 
the mass-mixing parameters. The result of such an exercise is shown in Fig. [6] which shows 
the {2a and 3a) ranges for values of f s > 0.1, where our Ly-a rescaling procedure furnishes 
robust results. 13 In the first panel (where the mass is projected out), the lower bound on 



(16) 



11 The log-normal distribution is more suited to treat the large hadronic uncertainties. 

12 Note that the three cases considered above (in which the relic abundance is univocally fixed by the 
mass and mixing parameters) are formally recovered introducing as penalty factor a Dirac delta function 
centered respectively on the^erage, lowest or highest abundance. 

13 It has been recently shown 41 1 that gamma-ray line emission searches can improve current X-ray limits 
in the region of very large masses and small mixing angles, which is relevant for very low values of the 
relic abundance (/ s < 0.1). 
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the mixing angle is determined by X-ray observations, while the upper bound is fixed by 
Ly-a measurements. The situation is reversed in the second panel (here the mixing angle is 
projected out) where the range allowed for the neutrino mass is limited from below by Ly-a 
measurements, and the upper bound is provided by X-ray observations. 

Although a direct measurement of the relic abundance of sterile neutrinos is not feasible, 
an indirect evidence for a mixed (WDM+CDM) scenario may constitute a realistic possi- 
bility. In this case, one might ask which parameters can give rise to sterile neutrinos in 
the right amount required to make the WDM component. Figure [6] provides an immediate 
and quantitative answer to this question showing at a glance where the DW model should 
"live" . As expected the allowed ranges gradually reduce for increasing values of f s until they 
eventually shrink in correspondence of the maximum value allowed for the relic abundance. 

We close this section with a final cautionary remark. As already stressed, the upper limit 
that we found on the relic abundance is valid only under the assumption that sterile neutri- 
nos are produced solely through the simplest production mechanism involving non-resonant 
oscillations. It can be evaded considering exotic models such as resonant oscillations 
boosted by large primeval lepton asymmetries, generation mechanisms operative during 



or prior to [49[ oscillations or low reheating temperatures [50|, |51[ • Therefore, further inves- 
tigation of these possibilities is important to prevent a viable dark matter candidate from 
being dismissed prematurely. In this context our results could serve as a quantitative guide 
for exploring possible mechanisms which are expected to act in synergy with the production 
via non-resonant oscillations. 



VI. CONCLUSIONS 

In light of recent results which seem to disfavor keV sterile neutrinos as viable dark matter 
candidates, we have revisited the constraints attainable on the parameters which govern the 
Dodelson-Widrow (DW) production mechanism. Relaxing the usual hypothesis that sterile 
neutrinos must account for all the dark matter content (Q s = fW)> we have shown how the 
X-ray and Ly-a measurements can be reinterpreted in the subdominant case (Q s < 0^). 
In addition, we have shown how the current data provide us with a conservative upper 
bound on the fraction of sterile neutrinos produced via the DW mechanism, which is robust 
with respect to the large uncertainties affecting the theoretical estimates. More sensitive 
X-ray observations, more quantitative Ly-a analyses, and a reduction of the theoretical 
uncertainties, all will need to play a crucial role in order to improve our limits. 
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FIG. 1: Constraints (at the 3a level) on mass and mixing parameters obtained from the spec- 
tral analysis of the HEAO-1 data under the assumption that sterile neutrinos account for all the 
dark matter content (O s = fidm)- The bounds obtained including only the Milky Way (MW) 
contribution, only the extra-galactic (EG) contribution, and their sum are shown separately. 
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Power suppression at k = 2h/Mpc 
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FIG. 2: Fractional suppression of the one dimensional matter power spectrum for mixed (CDM 
+WDM) models with respect to the case of pure CDM. The values refer to the pivot scale k = 
2/iMpc" 1 . The region under the thick solid line is excluded at 2a level by Lyman-a measurements. 
See the text for details. 
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FIG. 3: Comparison of the constraints on mass-mixing parameters imposed by experimental data 
(X-ray and Ly-a) with the theoretical predictions obtained in 75f] for the Dodelson-Widrow mech- 
anism. The two panels confront the "usual" case f s = l (left panel) with a representative case in 
which sterile neutrinos constitute only a fraction f s = 0.2 of dark matter (right panel). 



16 



10 



10 



X-ray 



average 
production 




Ly-a 



2cr constraints 
for 

three different 

theoretical 
estimates 



10' 



10 



lowest 
production 





_ i 2 -11 - 1 -9 
10 10 10 10 10 



io- 12 io- 11 io- 10 io- 9 10" 8 10" 7 



sin 2 (2a?) 



sin 2 (2a?) 



FIG. 4: Joint constraints from X-ray and Ly-a data on the mass-mixing parameters obtained for 
three different theoretical estimates of the relic abundance of sterile neutrinos produced via the 
Dodelson-Widrow mechanism. In each panel the allowed region lies under the solid curve, the 
dotted lines represent curves of constant sterile neutrino fractional abundance, while the dashed 
segment indicates its maximum value allowed in each case. In the first panel the "average" the- 
oretical calculation obtained in [75|] is used. In the last two panels the theoretical abundance is 



taken equal to the two extreme estimates determined in 75] 
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FIG. 5: Constraints on mass and mixing parameters governing the Dodelson-Widrow mechanism 
obtained by the joint analysis of X-ray and Ly-a data with the inclusion of the theoretical uncer- 
tainties. See the text for details. 
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